Influence of the field-detector coupling strength on the dynamical Casimir effect 
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We consider the problem of photon creation from vacuum inside an ideal cavity with vibrating 
walls in the resonance case, taking into account the interaction between the resonant field mode and 
a detector modeled by a quantum harmonic oscillator. The frequency of wall vibrations is taken to 
be twice the cavity normal frequency, modified due to the coupling with the detector. The dynamical 
equations are solved with the aid of the multiple scales method. Analytical expressions are obtained 
for the photon mean numbers and their variances for the field and detector modes, which are 
supposed to be initially in the vacuum quantum states. We analyze different regimes of excitation, 
depending on the ratio of the modulation depth of the time-dependent cavity eigenfrequency to the 
coupling strength between the cavity mode and detector. We show that statistical properties of the 
detector quantum state (variances of the photon numbers, photon distribution function, and the 
degree of quadrature squeezing) can be quite different from that of the field mode. Besides, the 
mean number of quanta in the detector mode increases with some time delay, compared with the 
field mode. 

PACS numbers: 42.50.Pq, 42.50.Lc, 42.50.Ar 



I. INTRODUCTION 

The effect of photon creation from vacuum in a cavity 
with rapidly varying geometrical or material properties, 
called nowadays as Dynamical Casimir Effect (DCE), at- 
tracted attention of many researchers for a long time: 
see recent reviews ,1 3]. The first experimental results 
for the open strip-line waveguide with time-dependent 
boundary conditions, simulating a single oscillating ideal 
mirror in the one-dimensional space, were reported re- 
cently in Q. It is quite probable that new experiments 
in other geometries will be done soon, too. Therefore, 
the problem of the back action of different detectors on 
the rate of photon generation becomes actual. The very 
first studies were performed in [j| Q under the condition 
that the field-detector coupling is much stronger than 
the effective coupling between the selected resonance field 
mode in the cavity and moving boundary (due to strong 
limitations on the attainable amplitude of surface oscil- 
lations in real materials). 

However, since that time there were proposed several 
schemes on simulating the DCE by means of changing 
not the positions of boundaries but their material prop- 
erties 0, 0] ■ I n such cases the effective velocity of bound- 
aries can be increased by several orders of magnitude, 
therefore the situations where the modulation depth of 
the instantaneous cavity eigenfrequency is bigger than 
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the normalized field-detector coupling coefficient can be 
considered now quite realistic, as well. 

There are two simple models of detectors. One of them 
describes the detector as a two-level "atom" It can 
be applied for the experimental setup proposed in 
This case was extensively studied in different regimes in 
[^-[lij], and generalizations to the three- and multi-level 
"atoms" were considered in (l3j . 

Here we consider another simple model, where the de- 
tector is represented by a harmonic oscillator tuned in 
resonance with the selected field mode 0, @| • This model 
seems to be adequate to the case of the so called MIR 
experiment @, [HI , where the microwave quanta created 
via the DCE are supposed to be detected by means of a 
small antenna put inside the cavity. Since the inductive 
antenna (a wire loop) used in that experiment is a part 
of a LC-contour, it can be reasonably approximated as 
a harmonic oscillator. Therefore we study the quantum 
system described by the Hamiltonian (in dimensionlcss 
units; in particular, we assume h = 1) 

H = -\pl + Lxj 2 (t)x\ +pl + oj 2 x 2 2 - Aojqkx^}. (1) 

Here x\ and p\ are the quadrature component operators 
of the selected resonance field mode, whereas x^ and p2 
are the quadrature component operators of the detector 
(antenna). The interaction part of the Hamiltonian is 
assumed to be proportional to the product X1P2 [H, Q in 
view of the standard minimal coupling term — (e/c)pA, 
assuming that operator x\ is proportional to the vector 
potential of the field mode. The instantaneous cavity 
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cigcnfrequency uj(t) is chosen in the form 
w(t) =uj [1 + 27COs(nt)], 



(2) 



where ljq is the unperturbed mode frequency and Q is 
the modulation frequency. The real modulation depth 
7 is assumed to be small, as well as the real coupling 
coefficient n. The influence of all other cavity modes can 
be neglected if the spectrum of cavity eigenfrequencies is 
not equidistant [TBj . 

If 7 = 0, then one can diagonalize Hamiltonian (JTJ, 
introducing two normal modes, whose frequencies are 
given by the exact formula lu± = ujq\/1 ± 2k, so that 
for \k\ C 1 we have 



(3) 



We assume that the modulation frequency O is twice 
bigger than one of the two splitted normal frequencies. 
We choose O = 2uj^ and normalize the frequencies and 
dynamical variables in such a way that wq = 1. Our 
goal is to study the influence of the dimensionless ratio 
(3 = j /(2k) on the number of excitations ("photons") in 
the field mode fii = \(x\ + pf) — \ and in the detector 
rt2 = \(x\ +P2) — \ for arbitrary (small or big) values of 
this parameter. 

The paper is organized as follows. In Sec. |TT]we derive 
the main dynamical equations and obtain their approxi- 
mate analytical solutions in the weak coupling and weak 
modulation regimes, comparing these solutions with re- 
sults of numerical calculations. In Sec. IIIII we use these 
solutions to calculate the time-dependent photon mean 
numbers for the initial vacuum state in both the cavity 
and detector modes. In Sec. IIVI we calculate the vari- 
ances of the photon number operator, the photon dis- 
tribution functions and the degree of squeezing in each 
mode. The results arc discussed in Sec. IVl 



II. MAIN DYNAMICAL EQUATIONS AND 
THEIR APPROXIMATE SOLUTIONS 



Hamiltonian ([T]) is a special case of generic quadratic 
Hamiltonians (see Ref. [l6| for a detailed treatment), 
which can be written as 

H = ^qBq, q = (p, x) = (pi, fa, x i: x 2 ). (4) 

Here B(t) is a symmetrical 4x4 matrix, which can be 
splitted in four 2x2 blocks as 



B 



h b 2 
63 o 4 



bi = bi, b 4 



where the tilde (") means matrix transposition. In the 
case involved we have 



61 



1 
1 



62 = 




-2k 



, 04 = 



LU 2 (t) 

1 



For any quadratic homogeneous Hamiltonian, the 
Heisenberg operators q(£) are related to the initial op- 
erators q(0) by means of a linear transformation: 

q(t) = £(i)q(0), $„(t) = C m (t)q a (Q) (5) 

(the summation over repeated Greek indices is implied), 
where coefficients £ Ma (i) form the time-dependent sym- 
plectic matrix C(t) = ||,£^ a (i)||. It is more convenient, 



however, to use the inverse matrix A 
([5]) can be rewritten as 

q (t) = A(t)q, 



£-\ Then Eq. 



(G) 



showing that qo(t) is the operator integral of motion in 
the Schrodinger picture. Matrix A(t) satisfies the follow- 
ing differential matrix equation and initial condition: 



A = ASS, 






h 


-h 






A(0) = h, (7) 



where In means the N x N identity matrix and the 
elements of matrix S satisfy the commutation relation 
£ij = i[xi,pj] = —Sij. Moreover, the following symplcc- 
tic identities are immediate consequencies of Eq. ([7J : 



ASA = E, A" 1 = EA" 1 = -EAE. 
Writing matrices A(t) and C(t) in the block form, 



A(t) = 



Ai A2 
A3 A4 



£(t) = 



A4 — A2 
—A3 Ai 



(8) 



(9) 



we obtain from Eqs. ([7|) the equations (for b\ = I2) 

dXi/dt = A163 - A 2 , d\ 2 /dt = A164 - A 2 6 2 , (10) 

and similar equations for the blocks A3 and A4 . Excluding 
matrices A2 and A4 we arrive at identical second-order 
equations for matrices Ai and A3: 



rf 2 Ai, 3 _ dXi^ 
dt 2 



dt 



■R-i + Ai^R3 — 0, 



where 



R, 



2k 



-1 

1 



R.3 



636 



3 02 



u 2 (t)-(2n) 2 
1 



(11) 



(12) 



(13) 



The difference between \\(t) and \%(t) is in the initial 
conditions: 

Ai(0) = 7 2 , A 1 (0) = 6 3 , 
A 3 (0) = 0, A 3 (0) = -7 2 . 

We suppose that I7I, |k| <C 1. In such a case, Eq. (jTTJ) 
can be solved with a sufficient accuracy analytically with 
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the aid of the method of multiple scales |l7|, [l8| . For this 
purpose we introduce a formal small parameter e, writing 
7 = £70, K = sko [with 70, Ko ~ 0(1)] and assuming that 
the solutions depend on the set of scaled times, To = u>ot, 
T\ = eujQt, T-2 = £ 2 uiQt, . . ., which can be considered as 
independent variables. This means that the time deriva- 
tives d/dt and d 2 /dt 2 can be expressed as follows: 



d d d 2 , 

dt = 0T +£ m +oi£) > 



dt 2 



d 2 

dT 2 



2s 



dTidTo 



0(e 2 



Also, we assume that matrices Xk (t) can be written in 
the form 

Afc(f) = A M (T , T u ...) + eA fcl i(T 0) T u ...) + 0(e 2 ). (14) 

For small enough values of parameter e it is sufficient to 
calculate the first term of the above expansion, taking 
into account the dependence on times To and T\ only. 
Therefore we neglect the terms of the second order with 
respect to e in Eq. (fTTj) . replacing exact matrix (| 13[) by 
the approximate form 



1 +47oecos(2T - 2k Ti) 
1 



(15) 



Omitting the details of cumbersome calculations, we 
bring here the general structure of solutions only. It ap- 
pears that the elements of each matrix Xk (t) (where 
k = 1, 2, 3, 4 and i,j = 1,2) can be written as (remember 
that we put Uq = 1) 

} = y cos(a;_t){ff > cosh^) + GW sinh(r M )/^ 



+U ( k l3) cos(tv) - 

^ sin ( w -*){/fc 
+u% j) cos(r t/ ) 4 



v k j) S1 



in(7v)/V} 



cosh(r /J 



to) 
9 k 



sinh(r M )/// 



to) 
w£ sin 



Or*)/"}- 



(16) 



The meaning of symbols is as follows. The "fast time" 
t appears in two oscillating functions cos(u;_i) and 
sm(co-t) only. The amplitudes of these fast oscilla- 
tions are modulated by combinations of hyperbolic and 
trigonometric functions, which depend on the "dimen- 
sionlcss slow time" t = nt (which is proportional to 
Ti). Namely, hyperbolic functions depend on variable 
= /it, whereas trigonometric functions depend on 
77 = vt. In turn, coefficients 77, /1 and v depend on 
the ratio j3 = 7/ (2k): 





= ^2ri + 2(/3 2 - 


1), 


(17) 




= v/2r ? - 2(/3 2 - 


1), 


(18) 


vW) 


= - 1) - 


f 1. 


(19) 



The factors fi 2 and v 2 are both positive for any values of 
j3 and satisfy the identity \w = 2/3. All other constant 



coefficients in Eq. ([16]) , such as , , and so on, 
also depend on the ratio /?. Explicit forms of coefficients 
with k = 1 and fc = 3 are given in Appendix. Since 
matrix 63 is of the order of e, we can write, in view of 
Eq. (fT0|) . approximate formulas A24 ~ —dXi^/dt. Con- 
sequently, within the chosen accuracy (i.e., neglecting 
terms proportional to e in all the amplitude coefficients) 
the elements of matrices A2.4 can be obtained from the re- 
spective formulas (fT()|) for A13 by means of replacements 
cos(u!-t) — > sin(o;_t) and sin(u-t) — > — cos(cj_£). 

We have checked the accuracy of approximate analyti- 
cal solutions (|16[) , comparing them with exact numerical 
solutions of Eq. (fTTj) . It appears that the coincidence 
is quite good for e < 10 -2 . An example of comparing 
analytical and numerical solutions is given in Fig. [1] 
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Figure 1. Numerical (points) and analytical (solid line) solu- 
tions for the matrix element (t, r) in the case of 70 = 4.0, 
k = 1.0 with p = 2 and e = 0.01. 



III. THE PHOTON GENERATION FROM 
VACUUM STATE 

Knowing matrices A(t) or C{t) ([9]) one can calculate 
immediately the time evolution of the symmetric co- 
variance matrix M. = |jA^ MI y||, whose elements Ai^ = 
h (q^q^ + qVqV) — (qV) are the central second-order 
statistical moments. Namely, 



M{t) = C(t)M(0)C(t). 



(20) 



For the initial vacuum states of both modes matrix Ai(0) 
is proportional to the unity matrix: A4(0) = Conse- 
quently, in view of Eqs Q and (|20[) , the mean number of 
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quanta in the mth mode (m = 1,2) can be expressed in 
terms of diagonal elements of the sum of products of ma- 
trices Afc and Afc (since the first-order mean values (q M ) 
are equal to zero in the case under study): 

1 /v-^~ \ (mm) 1 
2(X AfcAfe ) "2" (21) 

k=l 

These mean numbers do not depend on the "fast time" 
t. Their explicit expressions are as follows: 

*™ = ^2 { b + 1 + (! + u «) /? 2 (»7 + 2 )] § " 
+ [,,-l + (l+u m )/3 2 (r?-/3 2 )]S 2 
-2/3 2 (C Al C I ,-l) + 2/?u m S M S„}, (22) 

where ui = 1 and 112 = — 1. We have introduced the 
short notation 

C M = cosh(r M ), § M = sinh(r^), (23) 

C v = cos(t„), S„ = sm(r v ). (24) 

Typical regimes of excitations of the two modes are illus- 
trated in Fig. [2J where we plot the photon mean numbers 
in the cavity mode and the detector as functions of the 
"dimcnsionless slow time" r for different values of pa- 
rameter j3. Higher values of j3 imply on a more intense 
photon generation for a given time. 



We see that for j3 = 1 the numbers of excitations 
in each mode grow exponentially after a short transient 
time, but the number of quanta in the detector mode 
is significantly smaller than in the cavity mode for any 
fixed value of r. For smaller values of /3 we observe some 
"beats" between the two modes. However, these "beats" 
exist for limited time intervals, since only the first line in 
formula (|2"2"j) is important for ^> 1. Thus we see that 
asymptotically the number of quanta in the field mode 
is always bigger than in the detector (since 1 + Ui =2 
while 1 + u-2 = 0). The ratio can be very big if (3 >• 1: 
in this case rii/n2 ~ 4/3 2 = (7/ft) 2 . On the other hand, 
ni/"2 = 3 if /3 = 1, and Si rts V12 for (3 -c 1 (in the 
asymptotical regime r M 3> 1). 

For < 1 we have rj w 1 — (3 2 /2, v « 2, and fi « j3, 
so that t m w and 77, w 2k7j. Therefore the numbers 
of excitations in the both modes coincide (provided the 
slow time variable r is not very small), with corrections 
of the order of (3: 

n { £ <<1] = + f u m S, y § M + 0(/3 2 ). (25) 

Moreover, the main term in (f2"5]l does not depend on the 
field-detector coupling coefficient in this case. Eq. (|2"5)) 
is in agreement with the results obtained in Q by using 
the method of slowly varying amplitudes (although the 
argument of the sinh function in [6j was twice bigger than 
here due to a misprint in the definition of parameter /1 
in that paper). 

For /? = 1 we have r\ = 1, fx = v = y/2, and n = v^r. 
In this case 



3 



2 




2 3 
dimensionless slow time 



Figure 2. The behavior of function Y(t) = ln[l + n m (t)] for 
the cavity (1) and detector (2) with /3 = 0.2 (s), /3 = 0.5 
(/) and ,3 = 1.0 (/i). Note the time delay in increasing the 
number of quanta in the detector for /3 = 0.5 and /3 — 1.0. 



;(/3=l) _ 



= -[Si + Cx(d - d) + i^SiCSi + Si)], (26) 



so that ni 1 



(/3=i) 



3n 



(/3=1) 



if ri > 1. 



If /? > 1, then 7? « /3 2 



(j, « 2/3, and ^ w 1, so that 



7< and t v ~ nt. Therefore 



,(/3»l) 
l l 

;(/8>l) 



1 

" Z/? 2 



-C„) 



2sr 



0(/T 



(27) 
(28) 



The number of photons in the field mode coincides with 
the result [1, @] obtained in the absence of any detector, 
up to corrections of the order of 0((3~ 2 ). The number 
of excitations in the detector is about (3 2 times smaller, 
and this leads to a significant time delay in appearance of 
excitations in the detector: while the photons in the field 
mode appear after time t\ ~ 7 -1 , the detector begins to 
"feel" their presence after the time t-x ~ t\ m(/3). 

This time delay exists also for very small times, since 
in the limit r — > Eq. (f2"2"|) can be simplified as follows: 



ni =4/3 2 r 2 + 0(r 4 ), 
h 2 =/?V+ 0(r 6 ). 



(29) 
(30) 
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IV. PHOTON FLUCTUATIONS 



It is interesting to know, besides the mean numbers of 



quanta, their variances a„ = 



the degree of 



squeezing of quadrature components and the photon dis- 
tribution functions (PDF) V(m) = (m\p\m), where \m) 
means the mth Fock state and p is the statistical operator 
describing the mixed (due to the interaction) state of the 
mode. For the initial vacuum states of the held mode 
and detector the time-dependent statistical operator is 
Gaussian. Therefore we can use the results obtained in 
Ref. [l9| for the most general Gaussian states. 



A. Fluctuations of the photon numbers 



initial state of the total system is chosen to be pure (vac- 
uum), it remains pure for any instant of time in the ab- 
sence of dissipation (assumed in this paper). Therefore 
V 1 (t) = V 2 (t). 

Simple formula for variances can be obtained if (3 <C 1 : 



(2C 2 + 1)S* + ^S^C^C +4u m C M ). (35) 



r (/3«i) 



Comparing it with (|25p , we conclude that quantum states 
of the both modes are non-thermal, although they be- 
come highly mixed with the course of time, since the pu- 
rity in this case diminishes as Tr (p 2 ) = [cosh(7i/2)] _1 . 

The limit (f3"3")l is achieved approximately for the held 
mode if /3 ^> 1: 



In our case the first-order mean values of the quadra- 
ture components are equal to zero, so that 



r (/3»l) 



where 



2ni 



Mx r , 



2fL 



1 

4' 



Ma 



Mr 



(31) 



(32) 



is the invariant uncertainty product (IUP) of the mth 
mode, which satisfies the Schrodinger-Robertson uncer- 
tainty relation (l6j T> m > 1/4 for any quantum state. 
Besides, it is easy to verify that T> < (n + 1/2) for the 
states with zero first-order mean values of the quadra- 
ture components. Therefore the variances of the photon 
number in the Gaussian states with zero first-order mean 
values of the quadrature components must obey the in- 
equality er n > n(n+l). The equality sign is attained 
for thermal (i.e., mixed) quantum states. On the other 
hand, for pure squeezed vacuum states with T> = 1/4 we 
have 



cr„ = 2n(n+ 1) . 



(33) 



The calculations lead us to the identical expressions 
for the functions T> m (t) for m = 1 and m = 2 



V, 



^s 2 



+ (3rf + n 2 f3 2 - 2) (§ 2 - S 2 ) - 2/3 V 2 C^C„S„ 
+ 4/3 3 (C M C„-l)(§ Al S„-/3 3 )}. (34) 



This result may seem surprising at first glance, since 
other statistical properties (such as the mean number of 
quanta or degree of squeezing studied in the next sub- 
section) in each mode are different. But this is the con- 
sequence of the property of IUP for the Gaussian states: 
for these states IUP determines the quantum purity ac- 
cording to the formula Tr (p 2 ) = (42?) -1 / 2 . On the other 
hand, it is known that for any pure bipartite quantum 
state the purities of each part are identical. Since the 



2§ 2 C 2 



0{p~ 



(36) 



The invariant uncertainty product in this case equals 



i 

4~/3~ 2 



1 - 2C M C„] [1 + 0(/3- V )] 



(37) 

Comparing Eqs. (f28|) and ([37]). we conclude that quan- 
tum states of the both modes become significantly mixed 
when the mean number of excitations in the detector ex- 
ceeds the unit value. For r,, > 1 we have V m ~ exp(2r At ), 
whereas n 2 n ~ cxp(4r M ). Consequently, if n m » 1, then 
a nm ps 2n 2 n for m = 1,2, like in the squeezed vacuum 
state, even if the states are highly mixed. 

On the other hand, the photon statistics in the two 
modes are quite different for t m < 1, especially in the 
short-time limit r — > 0. Frequently the photon statis- 
tics is characterized by the Mandel parameter Q = 
(cr n — h) /n. However, in the case concerned this pa- 
rameter is not very useful, since it is always greater 
than unity and grows unlimitcdly with time in the DCE 
regime. A more convenient parameter is the ratio Z m = 
cr nm / [n m (1 + n m )], since it varies between 1 and 2 for 
the Gaussian states with zero mean values of quadrature 
components. The evolution of functions Z\{r) and ^(t) 
for different values of parameter f3 is illustrated in Fig. 
[3] The difference is clearly seen in the limit r — > 0, when 
Eq. (f3"4"]) can be simplified as follows: 



1 



rf 2 



l + -(3/? 2 -4). 



0(7 



(38) 



Therefore the contribution of T> m can be neglected in Eq. 
{3TJ with m = 1, in view of Eq. flU), so that Z X {G) = 2. 
But T>2 has the same order of magnitude as U2 for r — > 0, 
according to Eqs. ([3D]) and ((35J). Therefore Z 2 (0) = 1. 
More precisely, 



Zi(r) 



^2 r 



Z2(r) = l + - 



36 
19 s 

4 T 

45 



+ ^V+0(r 6 ), 
6 



6 



for ft = 1 wc obtain 




Xm 1/3=1 



3 + sin(2ri) 
4(2 + u m ) ' 



(41) 



This function oscillates between the values 1/6 and 1/3 
for m = 1 (the field mode). But it is three times bigger 
for the detector mode. It is not difficult to calculate the 
minimal value of the coefficient Xz as function of slow 
time tv for arbitrary values of parameter ft: 



v (as) l ■ 
A2 \mtn 



1 



PW - 1) 

2[l+i?C8)]' 



(42) 



Moreover, the minimum of this function with respect 
to ft also can be found analytically. It is achieved for 
Pmin = (\/3-l)/\/2~ 0.52, being equal to (l+y/3)(y/3- 
v / 2)/2 ~ 0.43. Consequently, the detector mode cannot 
be strongly squeezed. 



dimensionless slow time 



Photon distribution functions 



Figure 3. The time behavior of Z m = a Um / [n m (1 + n m )] for 
the cavity (1) and detector (2) with P = 0.2 (s), ft = 0.5 (/) 
and = 1.0 (h). 



B. Squeezing 

The degree of squeezing is characterized by the invari- 
ant squeezing coefficient (ISC) (20l | (which is equivalent 
to the principal squeezing introduced in [2lj ) 



4X> 



2h + 1 + v/(2n+l) 2 



4X> 



(39) 



This is twice the minimal value of variance of any quadra- 
ture component taken over the period of fast oscillations. 
For r M 3> 1 we have the asymptotical formula 



X 



(as) 



V/n m 

= V 2 (V + /? 2 ) + (3r/ 2 - 2) C 2 , + (2 - n 2 ) S 2 U 
- 2(3 (if - 2(3 2 ) C„S V 

x {2? 7 2 [/3 2 (7 / + /3 2 )(l+u„ i ) + l + r;]} _1 . (40) 



In the special case of f3 <C 1 we obtain x 



{as) 



1/2, in 



agreement with 0, @ . A very high degree of squeezing 

can be obtained in the field mode if j3 1: then Xi*^ = 
l/(4/3 4 ). However, there is practically no squeezing in 

this limit case in the detector mode: x^^ = 1 ~ C, y S^//3. 
For intermediate values of parameter ft the ISC does not 
go asymptotically to some limit value, but it exhibits slow 
oscillations in time with the frequency uk. For example, 



The PDF of the Gaussian states was derived in [19j, 
[221 . [23| . For zero mean values of quadrature components 
x and p it can be expressed in terms of the Legcndrc 
polynomials Pk(x) as 



V{k) = 



2Y_ 



fe/2 



1" 



(fe+l)/2 



Pi 



42? - 1 



where 



y± = l + 4P±2(l + 2h). 



(43) 



(44) 



The behavior of PDF as function of k depends on the 
value of the argument of the Legendre polynomial. If 
this argument is close to zero, then strong oscillations of 
function V{k) are observed, since Legendre polynomials 
of zero argument turn into zero for odd values of k. Oth- 
erwise V(k) changes slowly and monotonously. This hap- 
pens for ft <C 1, as was shown in If ft ^> 1, then the 
argument of the Legendre polynomials for the first mode 
tends asymptotically to the small quantity l/(4ift 2 ), so 
that the PDF shows typical oscillations of the squeezed 
vacuum state. But for the second mode the argument 
of the Legendre polynomials goes asymptotically to the 
value l/(iy/S), which is not small. Therefore the PDF of 
the detector mode is monotonous function without oscil- 
lations (that agrees with the absence of squeezing in this 
mode), which can be well approximated (for 1 <C k ~ n) 
by the universal dependence derived in [24j 



P{k) 



exp[-(2fc+l)/(4n)] 
y/ivn(2k+l) 



(45) 



In particular, considering formally k as a continuous vari- 
able (i.e., using the Euler-MacLaurin summation for- 
mula), we obtain the following value for the probability 
of detecting any value k smaller than the mean value n: 



Prob(fc < n) = crf(V2/2) w 0.68, 



(46) 
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where erf(z) = (2/y/W) J* exp(— x 2 )dx is the error func- 
tion. The distributions V\{k) and V^ik) for ft = 1 and 
r = 2 are shown in Fig. |4] They are rather different, 
especially for small values of k. The photon mean num- 
bers are fii = 112.4 and n2 = 38.7. Nonetheless, the to- 
tal probabilities to measure the photon number smaller 
than the n m are practically the same for each mode: 
Probi(/c < n) = 0.683 and Prob 2 (fc < n) = 0.681, in 
full agreement with Eq. (|46|) . 



if losses in the cavity can be neglected. But we leave this 
problem for the further studies, as well as the influence of 
different initial states (e.g., thermal or squeezed) of the 
field mode and detector on the rate of generation of the 
Casimir photons. 
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Appendix A 

Here we give explicit expressions for non-zero coeffi- 
cients in the solutions (fT6"|) for elements of A^-matrices 
with k = 1 and k = 3. 



F (ll) 


= v + 


P 2 , 


G[ 12) 


= -r? + l-/3(/3 


-1), 


_p(22) 


= v- 


P 2 , 


G (21) 


= TI-l + 0(p + 


1), 




= v- 


P 2 , 




= -r?-l + /5(/3 


-1), 


urn) 


= 7? + 


P 2 , 


y(2D 


= 77+1 -p(p + 


1), 



A (12) = 1 H 


-P, 




= P+1 


-(/3 2 +77)(2/3 4 


-1), 


A (2i) = _j 


-P, 




= P(P + 


l) + l-r?, 




4 12) = -1 


~P, 




= (P 2 ~ 


r))(2p + l)-p 


-1, 


4 21) = n 


-P, 


(21) 


= ~P(P 


+ 1) - 1 - V, 





Figure 4. The photon distribution functions V m (k) for the 
cavity (o) and detector (o) in the case of /3 = 1.0 and r = 2. 



V. CONCLUDING REMARKS 

We have shown that statistical properties of the de- 
tector quantum state (variances of the photon numbers, 
photon distribution function, and the degree of quadra- 
ture squeezing) can be quite different from that of the 
field mode. This can be important for planning the ex- 
periments on DCE and analysing their results. In par- 
ticular, the discovered time delay in the increase of the 
mean number of quanta in the detector indicates that 
influence of losses in the detector (which are not taken 
into account in the present study) can be essential, even 





= 1 


-P, 




= l-/?+(/3 2 + 77)(2/3 - 






= P 


-1, 


G (22) 


= p(p - 1) + 1 - v , 






= P 


-1, 


v 3 


= ( v - p 2 )(2p - 1) + p - 






= 1 


-P, 


t/(22) 
^3 


= -p(p-l)-l-T), 





f (U) _ 
J3 — 


-7/ 


-P 2 , 


(12) 
93 


= 77-1 


+ /3(/3 4 


■1), 


,(22) _ 
J3 — 


P 2 


- v, 


(21) 
93 


= 1-7/ 


-P(P- 


1), 


(11) 
«3 = 


P 2 


- V, 


(12) 
v 3 


= 77 + 1 


-P(P 4 


'1), 


(22) 
U 3 = 


-v 


-P 2 , 


(21) 
v 3 


= p(p- 


l)-7? 


- 1. 



All other coefficients of matrices Ai and A3 are equal to 
zero. The cocficicnts of matrices A2 and A4 are given by 

formulas = —fii\ f^A = ^13 > an d similar ones 

for other lower case and capital symbols. 
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